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Abstract. We study the entanglement spectrum of Heisenberg spin ladders of 
arbitrary spin length S in the perturbative regime of strong rung coupling. For isotropic 
spin coupling the entanglement spectrum is, within first order perturbation theory, 
always proportional to the energy spectrum of the single chain with a proportionality 
factor being also independent of S. In particular, although the spin ladder possesses an 
excitation gap over its ground state for any spin length, the entanglement spectrum is 
gapless for half-integer S and gapful otherwise. A more complicated situation arises for 
anisotropic ladders of higher spin S > 1 since here even the unperturbed ground state 
has a nontrivial entanglement spectrum. Finally we discuss related issues in dimerized 
spin chains. 



1. Introduction 

Quantum entanglement has by now developed to a key ingredient and tool of quantum 
many body physics[U [2]. More recently, the notion of the entanglement spectrum has 
provided a novel conceptual input to the field leading to new insights in the properties 
of various systems [3]. These comprise quantum Hall monolayers at fractional filling 
[HI SI El El d El El HOI [HI H2l H3l HH UHl H6], quantum Hall bilayers at filling factor 
v = lpl], spin systems of one [HI [191 [2Ql [2H [22l [231 12H E3 |26] and two [23 |28l [291 ED] 
spatial dimensions, and topological insulators [311 E2]- Other topics recently covered 
include rotating Bose-Einstein condensates [33], coupled Tomonaga-Luttinger liquids 
[31], and systems of Bose-Hubbard|35j and complex paired superfluidspj]. 

In Ref. [23] Poilblanc reported the observation that chain-chain entanglement 
spectra in two-leg spin-1/2 ladders are remarkably similar to the energy spectrum of a 
single spin-1/2 Heisenberg chain. Furthermore, the effective inverse temperature fitted 
to the data was found to depend on the ratio of the leg to the rung couplings and 
vanishes in the limit of strong rung coupling. Subsequently, one of the present authors 
observed a similarly striking resemblance between the entanglement spectra of quantum 
Hall bilayers at v — 1 and the energy spectrum of a single physical layer at half filling. [T7] 
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In the case of isotropic spin-1/2 ladders the above observations were explained only 
shortly later by Peschel and Chung [23] and by the present authors [26] via first-order 
perturbation theory around the limit of strong rung coupling. Moreover, as shown 
in Ref. |26j, an anisotropic spin-1/2 ladder will in general lead to an entanglement 
Hamiltonian with renormalized anisotropy, and in second order perturbation theory 
further corrections to the entanglement Hamiltonian occur. 

In the present paper we extend the above results for isotropic Heisenberg ladders 
to the case of arbitrary spin length S. We find that the entanglement spectrum is still, 
within first order perturbation theory, proportional to the energy spectrum of the single 
chain. This is an remarkable result since important features of the latter spectra, such as 
presence or absence of an excitation gap over the ground state, depend crucially on the 
spin length. [371 EE1 EH] Furthermore, the proportionality factor (effective temperature) 
between the entanglement spectrum and the single-chain energy spectrum spectrum 
turns out to be independent of the spin length. 

Moreover, we also investigate S = 1 ladders with uniaxial anisotropy. Here we 
encounter the situation that already the unperturbed ground state is not fully entangled 
leading to a zero-order reduced density operator which is, differently from the S — 1/2 
case studied earlier [26], not proportional to the unit matrix. 

This paper is organized as follows. In section [2] we analyze the entanglement 
spectrum of isotropic Heisenberg ladders of arbitrary spin length within first and second 
order perturbation theory around the limit of strong rung coupling. In first order, the 
entanglement spectrum is found to be proportional to the energy spectrum of the single 
chain, and the proportionality factor turns out to be independent of S. Our analytical 
findings are compared with results of exact numerical diagonalizations of small systems. 
Section [3] deals with spin ladders with anisotropic couplings and exemplifies the case of 
spin length S — 1. Technically complicated details pertaining to sections [2] and [3j are 
deferred to the appendices. In section [4] we remark on related observations in dimerized 
spin chains of higher spin length. We close with a summary and an outlook in section 

El 



2. Isotropic Heisenberg ladder of arbitrary spin 

We consider the Hamiltonian "H = "Ho + "Hi of an isotropic spin ladder of arbitrary spin 
length with 

^0 = Jr X] S2i&2i+l i (1) 
i 

"Hi = J\ X] {p2iS2i+2 + SW-l ^i+l) (2) 

i 

describing the coupling along the rungs and legs, respectively. The sites on, say, the 
first (second) leg are denoted by even (odd) labels such that the i-th rung consists of 
sites 2i and 2i + 1. All spin- 5* operators are taken to be dimensionless such that J r , 
Ji have dimension of energy. The following considerations apply to finite systems with 
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periodic boundary conditions as well as to those of infinite size. Let us now treat Tii 
as a perturbation to H with antiferromagnetic coupling, J r > 0. It is expected in 
general that the ground state of two leg spin S ladders are gapped for finite J r > 
and J\. For S = 1 for example this has been verified numerically in Ref. [3D]. The 
perturbative expansion around the strong rung coupling limit is therefore an expansion 
into an extended phase. The unperturbed ground state reads 



|0> 



\ s i) 



where the singlet state on each rung is given by, using obvious notation, 

_]\5-m 



\ S i) 



S 

E 



V2S + 1 



\m) 2i \ - m) 



2i+l 



(3) 



(4) 



In order to compute the ground state in first-order perturbation theory in Hi, let us 
first consider the states 

S2i\Si) — — S2i+l\Si) (5) 

which transform under general S'[/(2)-rotations generated by S 2 i + S21+1 as triplets. 
Thus, as a special case of the Wigner-Eckart theorem, the states S^|s;) and S'fjs^) are 
proportional to the triplet components \tf ) and |tf), respectively, which are obtained 
readily by expliciting Eqs. @,((5| and normalizing the results, 



\t]) 



(-l) s - m ^S(S + 1) - m(m + 1) 



S{3 + 1)(23 + 1) 



\t°i) 



K 1 ) 



s 

E 



x \m + l)ai 

Vs(-i) s - m 



m 



Sy/S(S+l) (2S+1) 



m) 2 i+i 



m) 



2i+l 



(6) 
(7) 



5-1 

E 



(-l) s - m y/S(S + 1) - m(m + 1) 



x \m) 



2i\ 



S(S + 1)(2S 
(m + l)) 2 i+i 



1) 



Considering now two neighboring rungs, the state 

S2iS2i+2\Si)\Si + i) = S2i+lS2i+3\Si)\Si + i) 



(8) 



(9) 



transforms as a singlet under arbitrary rotations and, according to the above finding, 
must be composed of triplet states on each rung. In other words, this state is 
proportional to (cf. again Eq. Q) 

(10) 



^(l^)ltfi)-lOl*° + i> + K rl >|^ + i>) • 



v/3 
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Indeed, a direct calculation yields the following matrix elements, 

<^ 1 |(*g 1 |^ atf2 |a i >| Si+1 > = - ls(S + l), (11) 

$|<tU&&*al*>l**i> =\S(S+1). (12) 

Note that the norm HS^S^M^Is^lsi+i)!! 2 = (S(S + l)) 2 /3 equals the sum of the squares 
of the above three scalar products, in accordance with the previous arguments. 

With these results, the correction to the ground state within first order perturbation 
theory can be formulated as 

|i) = ^^ + i)E[(---K 1 >K r + 1 i>---) 

-(•••I^)I4i)---) + (---I^ 1 )I4i>---)] ( 13 ) 

where the dots denote singlet states on each rung not explicitly specified. The reduced 
density operator is obtained by tracing out one of the legs from = (|0) + |1))((0| + (1|) 
and is given within first order in J\ / J r by 

2Ji 

(2S + 1) L \^ J r 

with L being the number of rungs. Again within first order perturbation theory, this 
result can be formulated as 



PS=77^^(1-^E^ +1 ) (14) 



Z 

with Z = tr exp(— "Hent) being a partition function with respect to the entanglement 
Hamiltonian 

%ent = E SiSi+l ■ (16) 

Remarkably, the prefactor 

= ^, (17) 

which can also be viewed as a formal inverse temperature, is independent of the spin 
length S. Fig. [T] shows the full entanglement spectrum of a S = 1 spin ladder obtained 
via numerical diagonalization in a wide range of Ji/J r . As seen in the right panel, in 
the limit of small Ji/J r the entanglement spectrum compares excellently to the energy 
spectrum of a single Heisenberg chain, as predicted by perturbation theory. 

In summary, the entanglement spectrum within first order perturbation theory 
perfectly matches the energy spectrum of a single spin chain. This observation has 
already been made in Refs. [24, 26J for S = 1/2 and is generalized here to arbitrary spin 
length. The latter result might be considered as somewhat unexpected since important 
features of energy spectra of spin chains, such as presence or absence of an excitation 
gap over the ground state, depend crucially on the spin length. j37J SHI EH] On the other 
hand, the energy spectrum of the underlying spin ladder has, in the limit of strong rung 
coupling, a gap over its ground state for any spin length S. 




Figure 1. (Color online) The entanglement spectrum of a S = 1 spin ladder with 
L = 8 rungs as a function of J r f J\ obtained via numerical diagonalization. In the 
right panel we have linearly rescaled the data in order to compare it directly with the 
spectrum shown on the very right of a single spin S = 1 chain of eight sites labeled by 
the ^-component of the total spin. The energies have been multiplied by 2, according to 
Eq. ( 16 ). For practical reasons the different sectors of total spin and linear momentum 
of the entanglement spectrum are not distinguished in the plot, but these quantum 
nubers perfectly match with those of the single chain. 



The contributions of second order in J\ / J r require lengthier calculations which are 
outlined in the appendix. The final result for the reduced density in up to second order 
is 

1 



o {2) 

Fred 



(2S + 1) L 



1 f2Ji 



i — -j- ^2 SiSi+i 



2 V Jr 



i 

I ((SiS i+1 ) 2 + - ^(S(S + l)fL 
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which can, within the same order, be reformulated as 



n {2) 

r red 



1 



exp ( —% 



/(2) 
'-ent 



(19) 



with Z = trexp(— % 

H-ent 



(2)- 



ent/ 



and 



2J/ 



i+l 



-5(5 , +i)x;^3 

5 ~r 



7 ^ (('S'i'S'i+l) (^+1^+2) + ('S't+l'Si+2) (<5i$s+l)) 

i £ ((&S+i) a + + + l)) 2 L 



For the case 5 = 1/2 the second-order entanglement Hamiltonian simplifies to 

1 rJr 2 



-1/(2) 
'T-ent 



2Ji 



^ S>iS>i+l — 2 J^J (^^+2 — SiS{ + i 



■ (20) 



(21) 



as already given in Ref. [26]. We note that the constant term in the Hamiltonian (20) 
is strictly speaking arbitrary since it is always canceled against contributions to the 
partition function. In Eq. (20) this constant was adjusted such that tr(%g n {) = 0. In 
Fig. [2] we compare the numerically obtained entanglement spectrum of an S = 1 spin 
ladder with L = 8 rungs with the energy spectrum of the second-order entanglement 
Hamiltonian (20). The plot shows nice agreement even if the coupling strength along 
the legs approaches closely the rung coupling. 



3. Ladders with anisotropic couplings 

The case of uniaxially anisotropic S = 1/2 spin ladders was already treated within 
first-order perturbation theory in the couplings along the legs in Refs. [2H 26J. As an 
important feature, the unperturbed ground state on each rung remains, for not too large 
anisotropy, the fully entangled singlet state. This is no longer the case for larger spins 
since the anisotropic contribution to the Hamiltonian couples the singlet to all higher 
multiplets with even total spin which will lead to a reduction of entanglement in the 
unperturbed ground state compared to the singlet state given in Eqs. @,Q. Indeed, 
it is easy to show by elementary means that for S > 1/2 the singlet Q is the only fully 
entangled state of two spins Su and S^+i which lies entirely in a single multiplet of the 
total spin. 

To give a specific example, let us consider the anisotropic rung Hamiltonian 

Ho = Jr^2 {^2i^2i+l + ^2i^2i+l + ^2i^2i+l) (22) 
i 

for S = 1. Here the ground state on rung i is given by 

\g t ) = a-\si) - ^j«+k°) (23) 

= 7s{{ a - ~ [l^l-^ + l- 1 ^ 1 ^] 
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Figure 2. (Color online) Cross symbols: numerical entanglement spectra obtained 
for an S — 1 spin ladder of L = 8 rungs. For each value of J r /Ji the levels have been 
rescaled to lie in the interval [0, 1]. Circles: rescaled energy spectrum of the analytical 



second-order entanglement Hamiltonian (20). Both data show good agreement even 
in the vicinity of J r j Ji ~ 1. 



a_ + 



V2a + ) 



|0>2*|0> 



2i+l 



(24) 



with 



\ 



1 / 3/2-(l-A)/ 6 

2 ' ' - /2 + A 2 /4 



(25) 



5 = 1 — A, and the quintuplet component \qf) being defined in Eq. (A. 6). This ground 
state leads to the unperturbed reduced density matrix 

\ 2 



o {0) 

r red 



Ay2 tt _ tt+ + ±clj s;s; | , (26) 



which is obviously not proportional to the unity operator but contains for A ^ 1 
anisotropic spin couplings of arbitrary range. In linear order in the deviation 5 from 
isotropy we have a + ~ Vz|5|/9, a_ ~ 1, and therefore 

1 



o {0) 

Fred 



+ o (5 2 ) 



(27) 
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or 



o {0) 

Fred 



with Z = trexp(— % 



(o). 



ent; 



and 



n 



(0) 
ent 



25 



£ ( s t s t 



(28) 



(29) 



Moreover, a more tedious calculation outlined in |Appendix B| yields the following result 
for the reduced density matrix in first order in Ji/J r and 5, 

2Ji 



Fred 



5 
3 



(30) 



Again, the presence of an only incompletely entangled unperturbed ground state leads 
to a considerably more complicated entanglement spectrum compared to isotropic spin 
ladders. Let us therefore return to the case of an isotropic Heisenberg Hamiltonian (fi]) 
along the rungs for again arbitrary spin length S. As a perturbation we introduce the 
bilinear Hamiltonian 



2i+l 



(31) 



describing an arbitrarily anisotropic coupling between next neighbors encoded, using 
sum convention, in the matrix A. Using the relations 



-3^i\si) = T\li8{8 + l)\ft 1 ) 



~$2i+l \ S i) 



\s(s + W) 



and 



tr 



2i+l 



tr 



2i+l 



WW 



W3/2 



S(S + 1)(2S+1 
v/3 



-°2i 5 



:5. 



2i 



(32) 
(33) 

(34) 
(35) 



'S(S + 1)(2S + 1) 

it is easy to see that the reduced density matrix in first order in Ji/J r is of the form 
(pb with 



-1/(1) 

/Lent 



2J; 

Jr 



(36) 
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Thus, the entanglement Hamiltonian in lowest nontrivial order perfectly reproduces the 
functional form of the perturbation along the legs. This result is easily generalized to the 
case of bilinear couplings of longer range, or spatially inhomogeneous couplings (with 



the matrix A depending on lattice indices). The above finding (36) can be seen as the 
consequence of conditions already formulated in Ref . [21] : (i) full entanglement of the 
unperturbed ground state, and (ii) the perturbation couples only to a single energy level 
(here the triplet) of the unperturbed Hamiltonian. 

4. Dimerized spin chains 

The above observation that spin ladders with isotropic couplings have perturbative 
entanglement spectra allowing for simple interpretations, can be generalized to spin 
systems of other geometry. As an example, let us consider a dimerized spin-5 1 chain (as 
opposed to a spin ladder) described by H = Ho + Hi with 

Ho = Jo X] SliSzi+l ' (37) 
i 

Hi — J\ y"l ■ (38) 

i 

For simplicity and definiteness we assume a even number L of lattice sites with periodic 
boundary conditions. Thus, the ground state \G ) (\Gi)) of Ho (Hi) consists of singlets 
coupling each even (odd) lattice site to the site with next higher label. Considering now 
Hi as a perturbation to Hq and tracing out all, say, odd lattice sites, one readily obtains 



a first-order reduced density matrix of the form ( 15 ) with 



?4nt — ~r ^2 S 2i S 2i+2 ■ (39) 



Here the missing factor of 2 compared to Eq. (16) just reflects that we consider here a 
single spin chain and not two-leg ladder. 

Moreover, it is also instructive to study a general linear combination of the above 
two singlet states, 



(< 



'cos(6/2)\G ) + 8in(0/2)e^|Gi>) (40) 

with two real parameters 9, <p. For S = 1/2 such states span the ground state space 
of the Majumdar-Ghosh spin chain. jUJ Tracing out again the odd-labeled sites from 
p — leads to 

Prod = (2g+ \ )L/2 (l + \sin9 (e-T+ + e*V)) (41) 

where the operator 




T = (g) \ \m) 2i+22i {m\ (42) 



translates all spin states by one lattice site of the reduced chain. Thus, the eigenstates 
\k) of p re d fulfill T\k) = e lk \k) where the wave number k is restricted to integer multiple 
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of An /L by the boundary conditions. In particular, the levels of the entanglement 
spectrum read as a function of k 

£(k) = - In (1 + sin 9 cos(y9 - k)) + ^ \n(2S + 1) . (43) 
5. Conclusions and outlook 

We have analyzed the entanglement spectrum of Heisenberg spin ladders of arbitrary 
spin length S in the perturbative regime of strong rung coupling. For isotropic spin 
coupling the the entanglement spectrum turns out, within first order perturbation 
theory, to be proportional to the energy spectrum of the single chain with the 
proportionality factor being independent of S. Additional corrections are obtained 
in second order perturbation theory which are also fully evaluated. Due to their 
proportionality, the energy spectrum of a single spin chain and the entanglement 
spectrum of the ladder share the property of being gapless for half-integer spin length 
and gapful otherwise. This is in contrast to the energy spectrum of the full ladder which 
is, for large rung coupling, always gapped for any S. 

A more complicated situation arises for anisotropic ladders of higher spin S > 1 
since here even the unperturbed ground state has a nontrivial entanglement spectrum. 
Finally we have discussed related issues in dimerized spin chains. 
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Appendix A. Second-order contribution to the ground state and the 
reduced density matrix of isotropic spin ladders 

Defining the two-rung singlet state 

4 } = ^ " IW> + l*. rl >l*J>) (A.i) 

the first-order contribution can be formulated as 

|l) = ^-^+l)E(-"^i---) • (A.2) 

The second-order correction to a non-degenerate state |0) is given by the general 
expression (using obvious notation) 

x (0|tti|0)(a|tti|0) 1. . ^ (OjjWWijO) 
|2) = £ ^ {Eo-E a )> - 2 £ {Eo _ Ea y 

The first term of the above r.h.s. does not contribute in the present case, and 
in the second term again only triplet excitations on the rungs occur which can be 
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evaluated using the matrix elements calculated previously. The last contribution is 
more complicated since here (for S > 1/2) also quintuplet excitations on a given rung 
come into play. The components of such a quintuplet on rung i read explicitly 

k 2 > =Q(S) E [(-l) s - m ^S(S+l)-m(m-l) 



m=-S+l 



+ 1) - m{m + l)\m + l) 2i \ - m + l) 2i+ i\ ,(A.4) 



kl) = Q(S) E [(-l) s - m ^S(S+l)-m(m + l) 

m=~S 

(2m + l)|m+l) 2i | -m) 2i+1 ] , (A.5) 



l??> =Q(S)\ - E [(-l) 5 - m (^+l)-3m 2 )|m) 2i |-m) 2m l (A.6) 

k" 1 ) = Q(S) E [(-i) s - m Vs(s + i)-m(m + i) 

m=-S 

• (2m + l)\m) 2i \ -m- l) 2i+1 ] , (A.7) 
l?r 2 > = Q(5) E [(-l) 5 - m v/^+l)-m(m-l) 

m=-S+l 

• v /,S'( 1 S'+l)-m(m+l)|m - l) 2i | - m - l) 2i+1 ] ,(A.8) 

where 

Q(5) = . V (A.9) 

yJS(S + 1)(2S + 1)(2S + 3)(2S - 1) 

The highest-weight state \qf) of the above multiplet is again readily obtained by, e.g., 
applying S 2i on \t]) and normalizing the result. Now in order to compute the second- 
order correction to the ground state one needs to evaluate expressions of the form 

(S 2i S 2i+2 + S 2i+1 S 2i+3 ) (• • • af] +l ■■■) . (A.10) 

If the rungs involved do not overlap, i.e. \i — j\ > 1, we just obtain another independent 
singlet <7j-*+i composed of triplets on rungs % and i + 1. In the remaining cases one finds, 
after somewhat tedious calculations, the following expansions, 

(S 2 i- 2 S 2i + S 2i -iS 2i+ iJ (• ■ ■ (TiJ + i ■ ■ •) 

V2 



, y S(S + 1)(2S + 3) (25 - 1) (• • • a^ i+1 ■ ■ ■) 

+ ^ + l)(- ffl ( Vr-) > (A-ll) 
{s 2i S 2i+2 + S 2i+ iS 2i+ ^ (• ■ ■ oj^+i ' ' ') 

- ^(2S + 3)(2S-l)(---a^ +1 -..) 

(•••*&. I-")" -45(5 + l)|0>. (A.12) 
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Here we have introduced two further types of singlet states, 



a 



(2) 



1 

7H 



a 



tjk 



HC%)> + IC 2 >kI>), 
kl^)^ 1 )- b£K> + 1^)14)) 



v/3 



v/3 



where 



14) 



It- } 



1 

To 

~3 



\t\M) 



3 

To 

/2 



i^kD+yfi^ 1 )!?!) 



IOk?> + 



:l*r x >W> 



l*i)l?7 ) ■ 



(A.13) 
(A. 14) 

(A.15) 

(A.16) 
(A.17) 
(A.18) 



are the components of a two-rung triplet state composed of a triplet on rung i and a 
quintuplet on rung j. Summing up, the second-order correction to the ground state can 
be written as 

JA 2 (S(S + 1)) 2 



|2> 



+ 



6 

(S(S + l)f 
6 

2V2 



-L\0) 

E (■ 

|i-il>l 



(i) 



r (l) 



' ' ' a j,j+l 



2(g(g + l)) ; 
3^ 



(7 



(1) 
i,i+2 



2V3V 



(i) 



(A.19) 



where the first line on the r.h.s results from the second term in Eq. (A. 3), and all other 
contributions stem from the second line there. With these results at hand, the second- 
order correction to the reduced density matrix is given by the sum of the expressions 

2 / \ 2 



trne g (|l)(l|) 



1 



J r J (2S + 1) L 



tr lleg (|0)(2| + |2)(0|) 



J^ 2 



~~ 5 E {(SiSi+i 



\ i 
1 

(2S + 1) L 

Si+lSi + 2 



(A.20) 



7~1 £>iSi + i 

i 

+ [Si+lS; 



i+2 



Si S. 
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1 / -> -> 2 -> -> 



13 



+ U(S + 1) £ 33+2 - ^(5(5 + 1)) 2 L 



(A.21; 



which lead to Eq. (18). The above considerations refer explicitly to the case S > 1/2 
due to the presence of the quintuplets which are absent for S = 1/2. In the second-order 
correction to the ground state rungs with such states occur in the second and the fourth 



line of Eq. (A. 19). The corresponding contributions to the reduced density matrix read 
J^ 2 



4/5 

J r J (25 + 1) L 

2 

~ 3 
2/3 



S(S + l)J2SiS i+2 , (A.22) 



Jr 2 



respectively, and vanish for S = 1/2 due to elementary identities of Pauli matrices. As 
a result, the S = 1/2 reduced density matrix in up to second order in Ji/ J r is given by 



o (2) 

r red 



1 

2^ 



J r V 2 V J r 



+ J ~ S>iS>i+l 



t+i 



16' 



(A.24) 



which was already found in Ref. [26] 



Appendix B. First-order perturbation theory for anisotropic ladders 

We now consider an Hamiltonian of the form (22) acting with amplitude Ji (cf. Eq. (|2 
as a perturbation along the legs. Its action on the each rung pair is given by 

( qX qX , ny ny I A Q z Q z 
\°2i°2i+2 ~r D 2i°2i+2 ~t~ ^ ^2? °2i+2 



5 oj+ N 



1 

3 



5 a + N 



(14)1^+1) + l^fei) 



+ 2 a +(i^)l^+i) + l^ki+i) 
Remarkably, no coupling to the excited rung eigenstates 

|ei) = + |^a-l^ ) 



(B.l) 



(B-2) 
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occurs, which are the other orthogonal states spanned by |sj) and |g°). Now the first- 
order correction to the ground state is readily obtained as 

,2- / «■ x2 



|1) 



Ji 2A 2 /3 

J T Eg + 2 

4 1/6 



A 



5 a + 



+ 



7 E 

Ji 3/4 
IE — 



1 



2a_ + 



5 a_ 



) E(---[K 1 >K r + 1 i) + K rl >M---) 



a 



E 



1*1 V2 



(B.3) 



where the dots denote states \gj) (cf. Eq. (23)) on all rungs j not explicitly specified, 
and 



E n = ~- 



1 1 - A 



A 2 

T 



(B.4) 



such that the ground state energy on each unperturbed rung is J r E g . With these results 
at hand it is straightforward to compute the reduced density matrix within first order in 
Ji/J r . For general anisotropy, however, the resulting expressions are rather cumbersome 
and difficult to interpret. We therefore concentrate here on the linear order in 5 = 1 — A, 



|1> 



J_i2_ 

'J t* 3 



+ 



Ji2 

J r- 3 



i-^)E(---K >l^ a >- 
1 + I s ) E (••• [\t])\t^) + iv 1 ) \t} +1 ) 



+ 5 



(B.5) 



where on each unspecified rung we have states of the form 

l*> = \s t ) - ^k°> + O (5 2 ) (B.6) 

which also build up the unperturbed ground state. Note that in linear order in 5 no 
contribution to |1) from quintuplets occurs. Finally, tracing out one of the legs leads 
one finds for the reduced density matrix 



o {1) 

Fred 



3 L 



25 



i-tE [Si % 



E 



SjS; 



i'Ji+l 



26 

T 



E 



QZ QZ 

j D i 



15 



i I qx qx I rara \ qz qz 

g \ i *+l » i+1 ) n i «+l 



75 



9 



+ 0{5 : 



(B.7) 



Using the identities S Z S Z S ± + S ± S Z S Z = S ± and (S z ) 3 = S z for S = 1 operators leads 
to the more symmetric form (30). 
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